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Abstract 

We use the geometric approach to the theory of Lie systems of dif- 
ferential equations in order to study dissipative Ermakov systems. We 
prove that there is a superposition rule for solutions of such equations. 
This fact enables us to express the general solution of a dissipative 
Milne-Pinney equation in terms of particular solutions of a system of 
second-order linear differential equations and a set of constants. 

PACS: 02.40.Yy 02.30.Hq 

1 Introduction 

An instance of nonlinear equations that has been receiving an increasing 
interest during the last years because of its ubiquity in physics and engineer- 
ing is the presently called Milne-Pinney equation [U [2] , 



k 



x = -u\t)x + — , (1) 

where /c is a real constant with values depending up on the field in which 
the equation is applied. Ermakov introduced this equation when looking 
for a first integral for the corresponding time-dependent harmonic oscillator 
[3j. In a very short paper, [2j, Pinney showed that the general solution 
of ([!]) can be written in terms of a pair of solutions of the corresponding 
harmonic oscillator and two constants. More recently it has been shown that 
the differential equation ([1]) admits a superposition rule involving any two 
particular solutions and two constants [1]. 



1 



Lately Haas [S] started studying how to find an approximate solution 
for the simplest damped Pinney equation, the one with a damping term 
linear in the velocity. Moreover other problems with a term with a quadratic 
dependence on the velocity have recently been studied O [71 [8]. 

The geometric theory of Lie systems [9]- [12] has been very efficient in 
dealing with equation ([T]) and the corresponding Ermakov system [21 IH [13] ■ 
The possibility of considering a larger class of systems that can be reduced 
to a Lie system has been shown recently. As a particular example we can 
study a generalization of the Pinney equation in which the constant k is 
replaced by a time-dependent function. So we can deal with an equation 
similar to ([T]) but including a term proportional to the velocity because such 
an equation can be reduced to one like ([ID in which the constant k is re- 
placed by a time- dependent term. The latter suggests revisiting the theory 
of Lie systems from this new perspective in order to deal with such dissipative 
Milne-Pinney equations. This is the aim of the paper, which is organized as 
follows: We study in Section 2 the possibility of removing the term linear in 
velocity by means of a time-dependent change of coordinates while the way 
of doing a similar thing by means of a time-reparametrization is analyzed 
in Section 3. The properties of Lie systems are used in Section 4 to estab- 
lish the corresponding nonlinear superposition rules for some time-dependent 
genealizations of Milne-Pinney equations. Finally in Section 5 the recently 
proposed theory of quasi-Lie systems is applied in the study of differential 
equations of such type. 

2 Damped harmonic oscillator with time-dependent 
angular frequency 

In this Section we show that we can remove terms proportional to ve- 
locity in a certain sort of second-order differential equation through a simple 
time-dependent transformation. As a main result we deal with a dissipative 
Milne-Pinney equation in order to remove its dissipative term in this way. 
Nevertheless this time-dependent change of variables is an ad-hoc method 
and we just explain here how it works. Later we will explain that quasi- 
Lie schemes explain this transformation [H]. We also obtain this result and 
more information about some dissipative Milne-Pinney equations through 
the theory of quasi-Lie systems without this ad-hoc assumption. 
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As a first example we consider an harmonic oscillator with a time-dependent 
angular frequency, uj{t), and a term proportional to the velocity with coeffi- 
cient 7(t), namely 

x + -f{t)x + uj^{t)x = 0. (2) 

Particular instances of this equation are the ones with a damping term with 
constant coefficient 70, 

X -|- 7o X + c<j^(t) X = , (3) 
or the interesting example, found in [15], given by the choice 



7(t) = -|log[/(t)], 

which describes a physical model with variable mass m{t) = [f(t)]^^. 

Consider a time- dependent transformation of the variable x of the form 

x = at)y (4) 

which leads to 

i = m y + m y, x = at) y + 2 at) y + cw y ■ 

This time-dependent change of variables transforms the equation (j2]) into 

H- 2 + .wf + .^w).^0. (5) 

Therefore we can eliminate the term proportional to velocity by using a 
function ({t) such that 

or, more explicitly, 

at) = Coexp(^-^l\{t')dt'^ . (6) 

In this case equation ([2]) reduces to the harmonic oscillator 

y + n'it)y = 0, (7) 
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with an angular frequency fl(t) given by 

m = ^Ht) - 2(^>! - i|> (8) 

because, in view of ({t) = -|7(t)C(t), we obtain 

C{t) 2 2 V 2 y 4 2 ' 

and then 

CW , .J{t)_ l\t) 7(t) 

CW ^^^CW 4 2 • 

Hence equation ([2]) is transformed into the time-dependent harmonic oscilla- 
tor (I7j) with the time- dependent angular frequency fl(t) given by (IHl). Note 
that equation ([71) is the reduced canonical form of ([2]) as it is indicated in 

We have found a time- dependent transformation which enables one to 
remove the term proportional to velocity in ([2]). We study some particular 
instances of this method. For instance, if 7(t) = 70, the transformation 

x = e-^'y, (9) 
transforms the differential equation ([2]) into 

u\t)-^y = 0. (10) 

So we can analyze the damped system ([3]) through the time-dependent har- 
monic oscillator ([7]) with a time-dependent angular frequency fl{t) given by 

m = ^\t) - 1 . 

This latter nonautonomous second-order differential equation has been con- 
sidered quite often both in the classical and in the quantum approach (see 
e.g. [HI [ini [20]) and we can deal with it by means of the theory of Lie 
systems [I3l [HI [21] . 

Had we started with a generalized Pinney equation with a time-dependent 
coupling k{t), 

X + -f{t) X + ij\t) X = , (11) 

x-^ 
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we would obtain, through the same time-dependent change of variables (jl]), 
the following equation 



at) J ' yat) "'cw "J ' cmy' 

Moreover, if ({t) is defined by ([6]), the preceding equation simplifies to the 
velocity-independent differential equation 

Here Q{t) is given by ([H]) and ( is the function ([6]). 

For the usually called damped Milne-Pinney equation, ■y(t) = 70, the 
transformation ([9]) reduces the differential equation ([3]) to 

Notice that in the particular case of equation ( fTTl) with 7(t) = 70 and 
k{t) = ko e"^"^"*, we recover a Milne-Pinney equation like ([I]) with ko and 
Q{t) instead of k and uj{t). This is exactly the example considered in [22] as 
associated with the Caldirola-Kanai model. 

In the more general case of 7(t) being an arbitrary function, if ({t) is 
chosen to be such that the term proportional to the velocity vanishes, i.e. 
k{t) = kQ exp(— 2 'j(t') dt'), we also recover a Milne-Pinney equation like 
([T]) with coefficient k^. 

To sum, for a generic function k{t) we can remove the term propor- 
tional to the velocity, but the coefficient of the nonlinear term becomes time- 
dependent. 



3 Time-reparametrization of some second-order 
differential equations. 

In the preceding section we studied a particular time-dependent change 
of variables which allows us to eliminate the term proportional to velocities 
in a certain kind of second-order differential equation. Now we show that 
this can also be done by means of a time-reparametrization. 
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Given the second-order differential equation, 

X = f{x,x,t), 

the time-reparametrization given by a new parameter s such that 

Tt ^ ^^^^ 

where a{t) has a constant sign, for instance a{t) is positive, which defines 
a good reparametrization allowing us to express s as a function of t and, 
conversely, i as a function of s, leads to 

X — — — — -— X <^=^ X — a[t) X 
ds a{t) ^ ' 

and therefore 

X = a{t) x! ^Q?{t)x" . 

So we get 

x" = (/(x, a{t{s)) x\ t{s)) - a{t{s)) x') . 

In the particular instance of a second-order differential equation of the 
type 

c{t) 



X — ait) X -\- hit) X -\- 

x"^ 



the transformed equation is 



I.e. 



OL{t{s)) \ >> a{t{s))J a'^{t{s)) a^{t{s))x^ 
Note that, if the function a{t) is chosen to be given by 



a{t) — ctocxp a{t') dt'^ , 



the term containing the new velocity disappears and an equation of Milne- 
Pinney type is obtained. Nevertheless the coefficient of the nonlinear term 
is s-dependent 

Ktjs)) citjs)) 

So, we have obtained a reduction process from the dissipative Milne-Pinney 
equation into an equation of Milne-Pinney type in a new way. 



4 Lie systems of second-order differential equa- 
tions. 

We remark that a system of second-order differential equations 

x' = f{x,x,t), z = l,...,n, 

is associated with a system of first-order differential equations with a double 
number of variables, namely 

— Ij'^ 

= f\x,v,t), ^ = 1,---,^- (14) 

Even more generally the system of n second-order differential equations cor- 
responding to the system of 2n first-order differential equations 

x' = a{t)v\ ._ , s 

V' = f{x,v,t), ^-^^■■■^n, U^J 

is related to 

« • 

x^ = —x^ + a fUx,x/a,t), i = l,...,n. 
a 

This type of equations may appear as Schrodinger equations of position- 
dependent mass with the von Roos prescription 

As particular instances, the damped time-dependent harmonic oscilla- 
tor described by ([3]) is associated with the system of ffist-order differential 
equations 

V = -e^ot^\t)x, ^^^^ 



and the damped generalized Milne-Pinney equation (1111) is related to the 
system 



X = e '^°*'v 



e... + • (17) 



Recall that the standard Milne-Pinney equation ([T]) is linked to the fol- 
lowing system of first-order differential equations 



X = V, 



uj\t)x + -^. 



k (18) 



X 



Hence we can study the second-order differential equation ([T]) through the 
system f|T8l) . The solutions of f|T8|) are integral curves for the time-dependent 
vector field 

d ( o. . k \ d 



A-W=„_+^-.^W. + -J_ (19) 

which can be written ClS db linear combination 

X(t) = X2 -cj2(t)Xi. (20) 

of the vector fields 

^1 = ^^' ^2 = ^^ + 4#- (21) 

ov ox X'^ ov 

These vector fields close on a real Lie algebra isomorphic to sl(2, M) with the 
vector field 

A-3 = if.#-.4V (22) 



2 \ dx dv ^ 

because the vector fields Xa satisfy the commutation relations 

[Xi, X2] = 2X3, [X3, XJ = -X2, [X3, Xi] = Xi . 

Therefore the second-order differential equation ([T]) is a SODE Lie system 
in the sense of fi3\ [2T] and we can use the theory of Lie systems in order to 
study its properties. 

More generally each time-dependent vector field of the form 

X(t) = Pit) Xi + a{t) X2 , (23) 
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where a and (3 are arbitrary time- dependent functions, is a Lie system of the 
same type as the one studied in [131 [21] . Its integral curves are the solutions 
of the system 

X = a{t) V, 

k (24) 

i) = l3(t) X + a(t) —. 

X"^ 

This system of differential equations is a generalization of the case given in 
[T3l [21] and is related to the second-order differential equation, 



X - - a{t)(3{t) X - a^{t)X 



a{ty 

because the derivative of the first equation is 

X = a(t) V + a(t)v, 



0, 



x-^ 



(25) 



and, in view of 



we obtain 



If a is positive and we use the notations a = e and f3 = —qe , equation 
([25]) becomes the more general differential equation studied in |25j . 

0. (26) 



X + I' X + qx ^ — 



The latter equation can be studied by means of the theory of Lie systems 
and a superposition rule for its general solution can be found. In order to 
get such a rule we have to consider two copies of the Lie system. 



X 



-F 

e V 



F 1 

V = —qe X 



(27) 



which corresponds to the second-order differential equation 

x + Fx + qx = 0, (28) 
Thus we obtain the system of first-order 



together with one copy of f[T8|) . 
differential equations 



X 

y 

z 



—F 



Vx 



k 

-ge^ X + e ^ 



-qe^y, 
p 

-qe 



x-' 



which corresponds to the differential equations of the integral curves for the 
time-dependent vector field 

F f d d d k d \ p f d d d 



\ ox Oy OZ X-^ OVx J \ OVx OVy OVz 

This vectorfield can be expressed as X = e"^ N2 — qe^Ni, where Ni and A^2 
are 

d d d d d did 

^^ = y^ — — ^2 = ^'x^ + + ^^^ + 

OVy OVx OVz OX ■ Oy OZ x'^ Ov^ 

These vector fields generate a 3- dimensional real Lie algebra with the vector 
field iVs given by 

^ 2\dx ^ dy ^ dz dv^ dvy dv^ 
In fact, they generate a Lie algebra isomorphic to s[(2,M) because 
[iVi, N2] = 2iV3, [iVg, iVi] = iVi, [iV2, iVa] = . 

The dimension of the distribution generated by these vector fields is three 
and the manifold of the Lie system has dimension six. There are three time- 
independent integrals of motion which turn out to be the Ermakov invariant 
Ji of the subsystem involving variables x and ?/, the Ermakov invariant I2 of 
the subsystem involving variables x and 2;, and the Wronskian W of the sub- 
system involving variables y and z. They define a foliation with 3-dimensional 
leaves. This foliation can be used to obtain a superposition rule. 

The Ermakov invariants are 

A = ^ {{yvx - xvyf + k ^ , ^2 = ^ {{xvz - zvxf + k 

where /i and I2 are non-negative constants for /c > and the Wronskian W 
is 

W = yvz- zvy . 

It is to be remarked that the relation of the first equation flTTj) allows us 
to rewrite the invariants Ji and I21 respectively, as: 

h = \ {e^^^'\yx - xyf + )') , h = \ [e'^^'\xz - zx f + k (^)') . 
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In particular, for F{t) = 70 1, which corresponds to the above mentioned 
example associated with the Caldirola-Kanai model, we recover the invariant 
given in (2.11b) of [22]. 

We can obtain an explicit expression of x in terms of y, z and the three 
first integrals Ii,l2,W 

X = ^ [l^y^ + hz^ ± ^/Ahh - kW^ yz) . 

We remark that is a constant fixed by the two independent particu- 
lar solutions of the time-dependent harmonic oscillator Xi{t) and X2it), i.e. 
y{t) = Xi(t) and z(t) = X2{t), and only Ji and I2 play the role of constants 
in this superposition rule for the Milne-Pinney equation. This is not a sur- 
prising fact because the Milne-Pinney equation is a second-order differential 
equation. Note also that the values of Ji and I2 are nonnegative constants, 
but should be chosen such that x(0) be real, i.e. 4J1J2 > kW'^. 



5 Ermakov systems that are not of Lie sys- 
tem type 

In this section we show that there exist dissipative Milne-Pinney equa- 
tions that can be transformed into simple Milne-Pinney equations by means 
of the methods developed for quasi-Lie systems. Furthermore we can deal 
with the so-obtained Ermakov systems as in [21] in order to find many of 
their properties, i.e. first integrals of motion or superposition rules. Next we 
can use these results to obtain properties of the initial system of differential 
equations by inverting the time-dependent change of variable. 
Consider the family of differential equations 

X = a(t)x + b(t)x + c(t)i- . (29) 

x-^ 

We are mainly interested in the case c{t) 7^ and we can assume that 
c{t) has a constant sign for the set of values of t we are considering. The case 
in which c{t) is identically zero corresponds to the harmonic oscillator with 
time-dependent frequency and a dissipative term. 

Usually we associate with such a second-order differential equation a sys- 
tem of first-order differential equations by introducing a new variable, v, and 
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relating fl29p to the system of first-order differential equations 



= a{t)v + b{t)x + c{t)i-. ^^^^ 

This system describes the integral curves for the time-dependent vector field 

X{t) = a{t)Xi + h{t)X2 + c(t)X3 + X4 , 
where the vector fields Xi, . . . , X4, are given by 

d d 1 (9 d 

Xi = v—, X2 = x— , X3 = — — , Xi = v—. (31) 
ov ov x'^ ov ox 

Consider also the vector field 

The set of these five vector fields is a basis for a M-linear space V. However, 
they do not close on a Lie algebra because the commutator [X3, X4] is not in 
V. Moreover, it can be checked that there is no finite-dimensional real Lie 
algebra V containing V. Therefore the differential equation (1301) cannot be 
considered as a Lie system. Nevertheless we can deal with this differential 
equation through a quasi-Lie scheme. 

The two-dimensional linear subspace, W G V, generated by the vector 
fields 

Fi = Xi = V—, F2 = X2 = X— , (33) 

ov ov 

is a Lie algebra because these vector fields satisfy the commutation relation 

[Y^,Y2] = -Y2. (34) 

What is more, as 

[^1,-^3] = —-^3, [li,X4] = X4, [Yi,X5] = 0, , s 

[F2, X3] = 0, [F2, X4] = X5 - Xi, [F2, X5] = -X2, ^''''^ 

the linear space V is invariant under the action of Lie algebra W on V , i.e. 
[PV, V] C V . Therefore we have found a quasi-Lie scheme to deal with the 
differential equation flSU]) . 
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The corresponding set of time-dependent diffeomorphisms of TM related 
to the flows of time-dependent vector fields in W, ai(t)Yi + a2{t)Y2, is given 
by 

X = x', 

V = a{t)v' + P{t)x' 
with a{t) 7^ 0. The inverse transformation is 



(36) 



X, 



m ^ 1 

-X H —V. 



(37) 



a{t) a{t) 



These time-dependent diffeomorphisms transform the system (13011 into a 
new one in which the time-dependent vector field determining the dynamics 
can be written as a linear combination of the fields of V at each time 



X'{t) = a'{t)Xi + h'{t)X2 + c'(t)X3 + d'{t)Xi + e'{t)X^ 
More explicitly the new coefficients are 

a{t) 



(38) 



a'{t) = a{t) - f3{t) - 
bit) 



b'it) 
c\t) 



a{t) 

a{t)' 
d'it) = a{t), 



+ a{t) 



aity 

m (3\t) Pit) 



ait) ait) ait) 



(39) 



[e'it)=Pit), 

and the integral curves for (1381) are solutions of the system 



dx' 

m 

dfv[_ 
lit 



(3it)x' + ait)v' 

m ( m 



+ 



ait) \pit) 

citl^ 

ait) 



+ ait) - Pit) 



Pit) 



x' + ait) - Pit) 



ait) 
a^t) 



(40) 
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However, notice that only, if (3(t) = 0, is this system associated with a 
second-order differential equation, more specifically with 



The Ermakov systems studied in [21j are in the family of differential 
equations fl30l) . Hence it is natural to look for sufficient conditions to be able 
to transform a given system of ( l30l) into one of these Ermakov systems of the 
form 

X = f{t)v, 

V = -u;\t)x + f{t)-^, ^^^^ 

where fc is a constant, which corresponds to the second-order differential 
equation 

X = ill X + m (^-uj'{t)x + . (42) 

Next we compare fHOll with (HTl) to transform equation (l30l) into one related to 
Ermakov systems. As a result we notice that a = f and the time-dependent 
coefficients a{t) and c{t) must be such that 

, / N c(t) 
ka{t) = --4, 

i.e. the sign of k must coincide with that of c{t) and 



This expression provides a sufficient condition in order to be able to trans- 
form a differential equation corresponding to the system fl5U]) into one of the 
form of (142|) . Taking the time-derivative on the first condition of (H3|) we 
obtain 2 k a a = c. Dividing by and using the second condition in (143!) we 
get 
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Thus the resulting transformation is determined by 



«w = y^, m = o, (45) 

for a certain constant k. 

Next we point out some differential equations ( 129|) appearing in the lit- 
erature that can be related by means of the method developed here with a 
particular Lie system: the Milne-Pinney equation [21] . 

As a first example we analyze the Chini differential equation [26] 

This equation is associated with the system of first-order differential equa- 
tions 

^ X = V, 

Pit) qit) ^ 1 k (47) 
X + 



2pit) pit) pit) x3 ■ 

This system is a particular instance of fl5U]) for the following choice of time- 
dependent coefficients 

ait) = -P^, bit) = -44> c(t) = . (48) 
Mi) P{t) Pit) ^ ' 

In this case, as c/(2c) = — p/(2p) = a, these coefficients satisfy the reducibil- 
ity condition (jUj) and we can transform this system into a Lie one through 
the transformation (136!) determined by the coefficients 

ait) = -i= , Pit) = , (49) 

i.e. by means of the time-dependent change of variables 

X = x' 



I 

V . 



(50) 



So equation fH7|) becomes 



dx' 




Ht 


Vpit) 


dv' 


qit) 


It 


Vpit) 



x' + 



(51) 



vW) 



x'^ 
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This system describes the integral curves for the time-dependent vector field 



X{t) 



1 



d 



V 



Introducing the time-reparametrization 

dt' 



'—+[-q{t)x' + — ] — 



X'- 



d 



dv' 



r{t) 



the system fl^T]) reduces to the Ermakov system studied in 



(52) 



dx' 

dr, 

dv' 



-q{t{T))x' + 



k 



(53) 



„/3 



dr ----- 
which corresponds to the second-order differential equation 



X' 



(54) 



Now we can use the Ermakov invariant, superposition rules etc. for (153|) . 
Next, by inverting the time-dependent change of variables used for the Chini 
equation, we obtain some kind of time-dependent superposition rule for the 
solutions of (H6i) . 

Walter [27] developed another interesting example: 



.. ^ Pit) ■ ^ Qit) 
X H — —X H — —X 



k 



pit) pit) p'^{t)x^ 



(55) 



This equation is associated with the system of first-order differential equa- 
tions 

' X = f , 

1 k (56) 



V 



Pit) q{t) ^ 

-V —X + 



p{t) p{t) p'^{t) x^ 
This equation is a particular instance of fl30|) with the functions 



Pit) 



Pity 



cit) 



k 



p'it) 



(57) 
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In this case, as c/c = — 2p/p, these functions satisfy the reducibihty condition 
f lH|) and if we choose 

a{t) = ^ /5(t) = 0, (58) 
pit) 

then the transformation fl5Bl) transforms the equation flS^ into the Lie system 

fix' 1 



dt p{t) 
Introducing a new time function r by means of 

we obtain 

|| = -g(t(r)Mt(r))x' + A, (61) 

which is a standard Milne-Pinney equation that can be described again 
through the theory of such systems as in [21j. 

Finally we consider the example of Colegrave and Abdalla [28] : 

x-2^^x + p\t)x = p\t)^^, (62) 

with an associated system of first-order differential equations 



X = V, 



' = -^v-p {t)x+p (t)— , 
p[t) x-^ 

which is a particular instance of fl29p with the coefficients 

G(t) = ^, 6(t) = V(t), cit) = kp\t). (64) 

In this case, if we take j3{t) = and we choose as a a solution of the differ- 
ential equation, 

||-oW=0. (66) 
17 



as 

ait)=p\t), (66) 
we obtain that (136!) transforms the initial differential equation into 

^ = At)v', 

which can be studied through the formalism of Lie systems. On the intro- 
duction of a new time function 

T{t) = j' p\t') dt' (68) 

the last differential equation reduces to 



dr^ p^{t{T)) x'3 
and the usual Milne-Pinney equation is recovered. 



^' + (69) 



6 Conclusions and Outlook 

We have studied ' damped time-dependent angular frequency harmonic 
oscillators' and time-reparametrizations of second-order differential equations 
through ad hoc transformations. As a result we have found a method to 
remove terms proportional to the velocity in certain set of second-order dif- 
ferential equations. Afterwards we revisited the theory of SODE Lie systems 
and quasi-Lie schemes. These theories allow us to study many (systems of) 
second-order differential equations. They provide methods to obtain time- 
dependent superposition rules, integrals of motion, solutions etc. In this 
paper we have applied these theories to some particular dissipative Ermakov 
systems to recover known properties from a new point of view. 

The theory of quasi-Lie schemes can be used to deal with many systems 
of differential equations. We expect to analyse new differential equations in 
forthcoming papers, recovering known properties and finding new ones. 



18 



Acknowledgements 



The authors are grateful to Profs. J.M.F. Bassalo, F. Haas and P.G.L. Leach 
for critical lecture of the manuscript and valuable comments. Partial finan- 
cial support by research projects MTM2006-10531 and E24/1 (DGA) are 
acknowledged. JdL also acknowledge a F.P.U. grant from the Minister io de 
Educacion y Ciencia. 

References 

[1] W.E. Milne, The numerical determination of characteristic num- 
bers, Phys. Rev. 35, 863-67 (1930). 

[2] E. Pinney, The nonlinear differential equation y" +p{x)y' -\-cy~^ — 
0, Proc. A.M.S. 1, 681 (1950). 

[3] V.P. Ermakov, Second-order differential equations. Conditions of 
complete integrability, Univ. Isz. Kiev Series 111 9, 1-25 (1880) 
(translation by A.O. Harin). See: Appl. Anal. Discrete Math. 2, 
123-48 (2008). 

[4] J.F. Carinena and J. de Lucas, A nonlinear superposition rule for 
solutions of Milne-Pinney equation, Phys. Lett. A 372, 5385-5389 
(2008). 

[5] F. Haas, Approximate solution for a damped Pinney equation, 
math-ph/0712.4083. 

[6] Z.E. Musielak, Standard and non-standard Lagrangians for dis- 
sipative dynamical systems with variable coefficients, J. Phys. A: 
Math. Theor. 41, 055205 (2008). 

[7] Z.E. Musielak, D. Roy and L.D. Swift, Method to derive La- 
grangian and Hamiltonian for a nonlinear dynamical system with 
variable coefficients. Chaos, Solitons and Fractals 38, 894-902 
(2008). 

[8] Shang-Wu Qian, Bo- Wen Huang and Zhi-Yu Gu, Ermakov invari- 
ant and the general solution for a damped harmonic oscillator with 
a force quadratic in velocity, J. Phys. A 34, 5613-5617 (2001). 



19 



[9] S. Lie and G. Scheffers, Vorlesungen iiber continuierliche Gruppen 
mit geometrischen und andcrcn Anwendungen, Edited and revised 
by G. Scheffers, Teubner, Leipzig, 1893. 

[10] J.F. Carinena, J. Grabowski and G. Marmo, Lie-Scheffers sys- 
tems: a geometric approach, Bibliopolis, Napoli, 2000. 

[11] J.F. Cariiiena and A. Ramos, A new geometric approach to Lie 
systems and physical applications, Acta Appl. Matfi. 70, 43-69 
(2002). 

[12] J.F. Cariiiena, J. Grabowski and G. Marmo, Superposition rules, 
Lie theorem and partial differential equations. Rep. Matfi. Pfiys. 
60, 237-258 (2007). 

[13] J.F. Cariiiena, J. de Lucas and M.F. Ranada, Recent applications 
of the theory of Lie systems in Ermakov systems, SIGMA 4, 031, 
18 p. (2008). 

[14] J.F. Carinena, J. Grabowski and J. de Lucas, Quasi-Lie schemes: 
theory and applications, ArXiv: 0810.1160. 

[15] LA. Pedrosa, Canonical transformations and exact invariants for 
dissipative systems, J. Matfi. Pfiys. 28, 2662-2664 (1987). 

[16] R. Mifson, Liouville transformation and exactly solvable 
Schrddinger equation. Int. J. Tfieor. Pfiys. 37, 1735-1752 (1998). 

[17] L.M. Berkovich, Transformation of Sturm- Liouville differential 
equations, Funct. Anaf. Appf. 16, 190-192 (1982). 

[18] T.f . Brazier and P.G.L. Leach, Invariants for dissipative systems 
and Noether's theorem, Rev. Mex. Fis. 40, 378-385 (1994) 

[19] J.F. Carificna, J. dc Lucas and A. Ramos, A geometric approach 
to time evolution operators of Lie quantum systems. To appear in 
fnt. J. Theor. Phys. 48 (2009) 

[20] J.M. Ccrvero and J. Viiiarroei, On the quantum theory of the 
dam,ped harmonic oscillator, J. Phys. A: Math. Gen. 17 2963- 
2971 (1984). 



20 



[21] J.F. Carinena, J. de Lucas and M. F. Rafiada, Nonlinear super- 
positions and Ermakov systems. In: Differential Geometric Metfi- 
ods in Mechanics and Field Theory, pp. 15-33, F. Cantrijn, M. 
Crampin and B. Langerock eds., Academia Press, Gent, 2007. 

[22] A. Nassar, Ermakov and non-Ermakov systems in quantum dissi- 
pative models, J. Math. Phys. 27, 755-58 (1986). 

[23] M. Aktas and R. Sever, Effective mass Schrodinger equation 
for exactly solvable class of one- dimensional potentials, J. Math. 
Chem. 43, 92-100 (2008). 

[24] O. von Roos, Position- dependent effective mass in semiconductor 
theory, Phys. Rev. B 27, 7547-7552 (1983). 

[25] R. Redheffer, Steen's equation and its generalizations, Aequa- 
tiones Math. 58, 60-72 (1999). 

[26] M. Chini, SuU'equazione del 2° ordine lineare omogenea, Atti. Ac- 
cad. Sci. Torino 33, 737-745 (1898). 

[27] J. Walter, Bemerkungen zu dem Grenzpunktfallkriterium von N. 
Levinson, Math. Z. 105, 345-350 (1968). 

[28] R.K. Colegrave and M. S. Abdalla, Invariants for the time- 
dependent harmonic oscillators, J. Phys. A: Math. Gen. 12, 3805- 
3815 (1983). 



21 



